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On the Bethe Ansatz for the Jaynes-Cummings-Gaudin model 
O. BabelorQ , D. Talalae\0 



Abstract. We investigate the quantum Jaynes-Cummings model - a particular case of the Gaudin 
model with one of the spins being infinite. Starting from the Bethe equations we derive Baxter's 
equation and from it a closed set of equations for the eigenvalues of the commuting Hamiltoni- 
ans. A scalar product in the separated variables representation is found for which the commuting 
Hamiltonians are Hermitian. In the semi classical limit the Bethe roots accumulate on very specific 
curves in the complex plane. We give the equation of these curves. They build up a system of cuts 
modeling the spectral curve as a two sheeted cover of the complex plane. Finally, we extend some 
of these results to the XXX Heisenberg spin chain. 



1 Introduction 



The Jaynes-Cummings-Gaudin model is defined by the Hamiltonian 

n— 1 n— 1 

j=0 j=0 

Here b, b^ is a quantum harmonic oscillator 

[6, 6 f ] = h 



and s?, Sj are quantum spin operators 



[st,sj]=2hs z j , [s z j ,sf] = ±hs 



j 



For the oscillator, we represent b, b' as 



b = h4-, tf = z (2) 

dz 



They act on the Bargman space Bb = <|/(z), entire function of z j \ f(z)\ 2 e ' R dzdz < oo j>. For 
the spin operators, we assume that s| acts on a spin Sj representation 



Sj\rrij) = hmj\mj) 



s^\rrij) = hJ Sj(sj + 1) — mj(m,j ± 1) \rrij ± 1), rrtj = —Sj, —Sj + 1, • • • , Sj — 1, Sj 
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where Sj is integer or half integer. 

The Jaynes-Cummings model is well known in condensed matter physics pQ El E] . It also appears 
in the book by M. Gaudin [5], but the connection between the two seems to be recent [I]. 



2 Bethe Ansatz 



In order to write the Bethe Ansatz, we introduce the Lax matrix 

M(A) B(X)\ 
L{A) -\C(X) -A(X)J 

where the operator valued matrix elements are defined as 

A(X) = - + > 

g 2 2 Z^X-ej 



3=0 



2b ^ a; 



9 U X -*> 



9/it n ~ l S + 

c(x) = ^ + £ s > 

3=0 



9 ~ A - e j 



It is simple to check the commutation relations 

[A(X),B(n)} = *(B(\)-BM) 
A — fx 

L4(A),C( M )] = -_5_(C(A)-C(m)) 

[B(X),C(f,)} = ™(A(\)-AM) 
A — /x 

Moreover one has [A(\),A(p)] = 0, [B(\),B(p)] = and [C(\),C(p)] = 0. Defining 

^Tr (L 2 (A)) = A\X) + ^(-B(A)C(A) + C(A)S(A)) 

We have [Tr L 2 (A), Tr L 2 (/i)] = so that TrL 2 (A) generates a family of commuting quantities. 
Expanding in A we get 

1™ /^9^xn 1 n9 4 tt 2 ^ Hj ^ h 2 Sj(sj + 1) . . 

-T r (£»(»)) = ? (2A - „)> + ? if„ + - 2 g ^- + g ( /_' e . )2 (3) 

The Hamiltonian eq.([T]) is given by 

3 
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To write the Bethe Ansatz we define the reference state which is the lowest weight vector: 

|0) = |0)®|-si)®---®|-s n ), 6|0>=0, sj\- Sj )=0 
This vector has the following important properties 

fl(A)|0>=0 

and 

A(A)|0)=a(A)|0), a{X) = 2 4-- 2 -Y,^ 

Moreover since [5(A), C(A)] = 2hA'{\) we also have 

B(A)C(A)|0) = 2ha'(X)\0) 

With all this we deduce 

iTrL 2 (A)|0) = (a 2 (A) + W(A))|0) 

Let us now define the vector 

fi(A*i,M2, • • • = C(Mi)C(M2) • • • C(hm)\0) 
It is not difficult to prove that (see e.g. [5]) 

-TrL 2 (A)J2(Ml>M2, • • -,Hm) = A(A,^i,/i 2 , • • • , hm)&{hu (Jq, • • • , MAf) 

A(A, Ml, ^2, • • • , Mm) = a 2 (A) + ha'(X) + 2^ ° (A j ~ ^ (4) 
provided the parameters /ij satisfy the set of Bethe equations. 

3 Riccati equation 

We now analyse the Bethe equations eqs.flSJ). We introduce the function 

. z Hi 

Proposition 1 The Bethe equations (0) imply the following Riccati equation on S(z) 

S'(z) + S 2 (z) + A ((2* - u)S(z) -2M)=J2 2sj S{z) l S{6j) (6) 
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Proof . The Bethe equations read 



"lit, to hsj ^ h 



g 2 g 2 ^ Hi- €j <rf. m - Hj 

we multiply by l/(z — Hi) to get 

2 Hi u 1 ^-^ hsj 1 ^-^ h 1 



a , - ii a ii - ii 



9 2 z- ^ g 2 z- ^ <-f Hi ~ e j z ~ Hi *rf. Hi~ Hj z - Hi 
We now sum over i. We have 

M 

y_H^ = yHi^ + ^_ = _ M + zS{z) 
^ z-m z- m z - m 

also 

■ ' 1 1 1^/1 1 \ S{z) - S(ei) 



J2— — - =— £ — + — 



Z — €j 



and finally 

AI 1 1 1 M 1/1 1 \ 1 ^ 



. 4 . in- Ho* -Hi 2 ; { f^,i i; /'./ v- /'- • /'// - f^jr; (: /''' (: - 



In equation ([6]) the S'(ej) appear as parameters. The Riccati equation itself determines them as we 
now see. Suppose first that Sj = 1/2. We let z — > e« into eq.© getting 



S'(e<) + 5 2 ( £i ) + ((2e, - a;)^) - 2m) = S'fe) + £ 



%r v / rrf e» - ej 

The remarkable thing is that S'(ei) cancel in this equation and we get a set of closed algebraic 
equations for the S(ej). 

^ + h (>* ~ ~ 2M ) = E S(ei) Z S(€j) » i = l,-,« W 



Suppose next that Sj = 1. We expand the Riccati equation around z = ef. 

(z - e,) = S'ia) + S 2 (e l ) + A(( 2ei - a;)^) - 2M) = 2S'(q) + 2 ^ 5(ei) " ^ 



lig e i ~ e j 



3+i 



(z-ei) 1 : S"(e t ) + 2S(e i )S'(e i ) + ^((2e i -uj)S f (e l ) + 2S(ei)) 



S"(6,)-2£ 



S{€j) - Sjej) S'(ej) 
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We see that in the second equation S"(ei) cancel. The first equation allows to compute S'(ej) and 
the second equation then gives a set of closed equations for the S(ei). The general mechanism is 
clear. For a spin s, we expand 

s , _ 2s S(z)-S{e) = y ^2s s (m) {e){z _ sm-l 

z — e z — ' ml 

m 

and we see that the coefficient of S^ 2s \e) vanishes in the term m = 2s. The equations coming from 
(z — e) m_1 for m = 1, • • ■ , 2s — 1 allow to compute S'(e), ■ ■ ■ , S ( - 2s ~ 1 \e) by solving at each stage 
a linear equation. Plugging into the equation for m = 2s , we obtain a closed equation of degree 
2s + 1 for 5(e). 

iWS(e)) =0 (8) 

Notice that if M < 2s, the system will truncate at level M because there always exists a relation 
of the form = P(S, S', ■ ■ ■ S^- 1 )). 

The S(ej) also determine the eigenvalues of the commuting Hamiltonians. Going back to eq.Q, 
we see that 

a(A) — a(/ij) 2 Hsj ( 1 1 



X- m g 2 ^X-Hi\X — ej [ii 



2 hsj 1 _ 2M t-^ hsjS(ej 



\ _ 7 „, _ 7 = ~7f~ ~ 2^ 



Hence 



A(A) = a 2 (A) + fia'(A) +2hh^-"£ j (9) 

Expanding a(A) we deduce the eigenvalues of the commuting Hamiltonians 

h n = HM - Sj + - , = -2-% - -^^ej - 2/i SjSiej) + 2 ^ — ^7- 
j 9 9 ^.e, e; 

The algebraic equations for the S(ej) are therefore the characteristic equations of the set of matrices 
Hj. The existence of such characteristic equations is a general phenomenon. In the Appendix we 
derive them for the Heisenberg XXX spin chain. 



4 Baxter equation 



We can linearize the Riccati equation eq.© by setting 



Obviously 

M 



Hz)=n(z-^) (11) 
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The linearized equation reads 



r(z) + j.a(zW(z) + - 1-^- + Y, z - €j J ^ = ^ 

Here, we should understand that the S(ej) are determined by the procedure explained in the 
previous section, for instance by eq.fjZj) for spins Sj = 1/2. For such values of the parameters, the 
equation has the following remarkable property. 



Proposition 2 (Mukhin, Tatasov, Varchenko |6j) For the special values of the parameters 
S(ej) coming from the Bethe equations, the solutions of eg. fflty have trivial monodromy. 



Proof . Strictly speaking the proof in [B] is valid for the finite-dimensional representations. We 
provide here a straightforward proof in our case. The proposition is clear for the solution ipi{z) 
defined by eq. ([TT|) since it is a polynomial. A second solution can be constructed as usual 



exp 



2 f y 

h 



a(t)dt - 2 log V>i (y) dy = ^i(z) 



f < j % 



dy 



The monodromy will be trivial if the pole at y = Hi has no residue preventing the apparition of 
logarithms. Expanding around Hit we have 



exp , -- 



2 r y 



a{t)dt- 2 log -01 (y) 



exp 



h 



(y - h%Y 
(y - im) 



Hi ~ Hj 



+ 0(y - m 



but the coefficient of the dangerous (y — Hi) term vanishes by virtue of the Bethe equations. 



Next we set 

iP(z) = exp (-1 [ Z a{y)dy ) Q(z) (13) 



n 



We obtain for Q(z) the equation 



h 2 Q"(z) - \a 2 (z) + ha'{z) + *™ - 2h 2 £ j Q{z) 

Comparing with eq.©, this is also 

h 2 Q"{z) - K{z)Q{z) (14) 
Hence, we have recovered Baxter's equation. Notice that 

■^{P-wz) -k I (z 2 ~LUz) M 

Qi ^i[^*^w^W- IM) <15> 
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5 Bethe eigenvectors and separated variables 



We recall the form of the Bethe eigenvectors. 

n(fj,i,fj, 2 , • • • , hm) = c(m) ■ ■ ■ c(hm 

Following Sklyanin, [7J, we introduce the set of zeros of C(A) 

2*nr=i(A-A,) 



C(A) 



9 n-=i(A-e,) 



The operators Aj commute among themselves. Inserting this expression for C(A) into the Bethe 
state and remembering eq. (jlip . we find 

If we now switch to a Schroedinger representation where the Aj are represented as multiplication 
operators, the eigenstate eq. (fT6|) is represented by a product of functions in one variable tp(\i). 
This means that we have separated the variables in the Schroedinger equation. 



Proposition 3 In the separated variables, the Hamiltonians read 

llfc( e j-Afc) s^Tik^Mi-tk) ( d 2 2 d M 



Hj ~ Ihrtb ~ e k ) \ n^(A. - A fe ) Ua| + h^dK h I (17) 



Proof . Let us introduce the set of commuting operators Hj diagonal in the Bethe states basis 
eq. ()16p . (here we assume completeness of Bethe Ansatz), and such that 

HjH^Xk) = 2s J S(e j )H^(X k ) 

k k 

These are essentially the same operators as in eq.([3]). Then eq. (|12p implies for each variable Aj 

j k k 



Since this formula holds for a basis of eigenvectors, we can "divide" by Ilfc^KAfc). Inverting the 
Cauchy matrix Bij = l/(Aj — €j), and taking care of the order of operators we obtain the Hamil- 
tonians Hj in terms of the separated variables 

, 1 d 2 2 , d M 

explicitly, they are given by eqs. (|17j) . These Hamiltonians are known to commute [5JE1QI]]. ■ 
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To be able to work in this representation, we need the scalar product. We set 

IM| 2 = JUdKdKWW p(x 1 ,x 2 ,---,x n )H\^(X l )\ 2 (18) 



where 

W = H(Xi - Xj 

and 

_ 1 n 3 -|e«-Aj 



The measure p(x\,X2, ■ ■ ■ ,x n ) is determined by requiring that the Hamiltonian Hj are Hermitian. 
Proposition 4 The Hamiltonians Hj are Hermitian with respect to the scalar product eq. U8\) if 

I \ f°° ,1 -V M+n-Y-(si+l/2) TT Jzsi+l&y/yXi) , , 

p(xi,x 2 ,---,Xn) = dye y y + ^ ^ ' > || (19) 

where J2 Si +i{x) is the Bessel function. For n = 1, the formula for p{x) can be simplified giving 

p(x) = +1 [e- x x M ~ 2s ] = e- x P M -2s(x) 
where Pm-2s( x ) is a Laguerre polynomial of degree M — 2s. 

Proof . We have to show that 

f Y[dX k dX k \iP(X k )\ 2 £ ( + | A • a(A J ) + ^ j B^\ W \ 2 p(x u ■■■,x n ) (20) 

k j \ 3 ■? / 

is real. Now l?^ 1 = A _1 Ajj where Ajj is the minor of the element Bji. It is clearly independent of 
Xj. Hence 



dXj 13 « \dXj dXj 



We have 



where we introduced 



3%; 



These variables satisfy -^.z k = Bj k z k so that A _1 ^-A = W~ 1 -^-W — Bj k and therefore 
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Next, we have 

—p(xx,---,x n ) = ^2B jk x k —p(xi,---,x n ) 

3 k 

d 2 d 2 
-^p(x lr --,x n ) = B jk B jlXkXl Qx Qx p( Xl , ■ ■ ■ , x n ) 

i k,l k i 

= yZ^ffc^fc^TP^i)*"'^) + 2 y] S ifc — " — gfc^ o p(xi,---,x n ) 
J ax% e k -ei dx k dx\ 



Putting everything together eq. (|20j) becomes 

J HdXtdX^X^lWfY^^ 1 I -J2B%x k V k + Y, B Jk°k + \ V * } P(zi,-' -,x n ) (21) 



k k 



where we have defined 

V k = x k d 2 k + 2{s k + l)d Xk V Q = Y^ x k d Xk + M + n + 1 

and 

O k = t ( e fc - 77) - 2 — - — {{x k d Xk + s k + l){xid x + St + 1) - s k si - 1 

The conditions on p(xi, • • • , x n ) are that the sum over j in eq. ([2TJ) should be equal to its complex 
conjugate. When we perform this sum, we first get 

Y J B^B ]k O kP = O iP 

j,k 

which is real and gives no condition. Next we have the identities 

3 

1 z k 



The conditions on • • • , cc n ) then read 



-{zi - Zi) [V iP + V 0P ] + k i _ k i [D iP - V k p\ = 

k+i ei €k 
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Finally we find the n conditions 



V iP -V kP = 0, k + i (22) 
V iP + V p = (23) 



Notice that eg. (|23H is independent of i if the conditions eq. (122p are satisfied. A solution of eg. ([22 
is 

oo n 

p{x l ,---,x n ) =^2c p n 

p=0 5i H \-qn=p 



~9« 



*<ft!(2*i + l + fc)! 



Then eq. (|23j) gives 
the solution of which is 

Hence we have found 



Cp+i + (M + n + p+ 1)C P = 



r„ ( M ^ - /; ) : 



^ 1 ,x 2 ,... ! x n ) = ^(-ir( M + n + np! J; II — . .... < 2 ') 



This is equivalent to eq. (|T9|) . 



This important formula should be further studied. In particular, for ip(z) being a Bethe state, one 
should be able to compute it exactly because we know that by Gaudin formula 

||n(/ii,/i 2 , • • • , /xm)|| 2 ^ det J 
where J is the Jacobian matrix of Bethe's equations. This is still very mysterious. 

6 Semi- Classical limit 



The exact formula relating Q(z) and ip(z), eq. (fT5j) . allows to study the properties of the solutions 
of Bethe roots fj,i in the quasi-classical limit h — > 0. Let us set 



y(z) = h- 



Q(z) 

Then Baxter's equation, eq. (|14[) . becomes 

hy'(z)+y 2 (z)=A(z) (25) 

where A(z) is defined in eq.([9]). This is just another form of eq.([6]). In the semi-classical limit 
eq. ([25p becomes the equation of the spectral curve of the model (in that limit hsj = O(h )): 

y\z) = A(z) 
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/,From eq, (|T3|) we deduce that 

y(z) = a{z) + V 

so that we expect in the semi-classical limit 

. z fj, t 

This is a remarkable formula. It gives us the distribution of Bethe roots /x, in the semi-classical 
limit, as we now show. Let y/ A(z) be represented as a meromorphic function in the cut z-plane. 
Let us put the cuts so that 

2z ijj 



\J A{z) = - -z + 0{z 1 ), |z|^oo 
9 9 



and (we neglect terms of order h which do not contribute in the leading h approximation). 



VMz)= -J_ + 0(1), Z^€j 



By Cauchy theorem, we have 



dz' y/Ajz 7 ) 



lc 2m z' — z 

where C is composed of a big circle Co at infinity, minus small circles Cj around z = ej, minus 
contours Ai around the cuts of yj A(z). Hence 



r—- r dz' TaR) ^ r dz' ^J) ^ r dz' y/Ag) 

V [Z> 2irr z' - z ^ r 2m z' - z ^] A .2m z'-z 



C 



But 



and 



Co 2nr z'-z g l g 2 



dz' yj A(z') f dz' —fr- s j fc s j 



so that we arrive at 



c . 2m z'-z J Cj 2m (z' - ej)(z' - z) z- ej 

dz' y/Ajk 1 ) 



y/XJzj - a{z) = -J2 I 

: J A 

and therefore we should identify 



A. 2m z' — z 



z — Mi J Ai 2Z7T z — z 

Comparing both members of this formula suggests that the Bethe roots \ii accumulate in the semi- 
classical limit on curves Ai along which the singularities of both side should match. To determine 
these curves we assume that the Bethe roots //j tend to a continuous function n(t) when H — >• ( 
t = hi and % = OifT 1 )). 



y^ = y h ^ f dt = f dJ^-)^- 

^z-m ^z-n(i) J z-n(t) J A \dfij z- 



/' 
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Here A = YL-^i- Hence, comparing with the semi-classical result, we conclude that the function 
/j,(t) should satisfy the differential equation 



dt yW) 

The boundary condition is that the integral curve fi(t) should start (and end !) at a branch point 
of the spectral curve y 2 = A(z). We stress that the function A(z) is completely determined by the 
Bethe equations themselves so that these equations "know" the Riemann surface. 

This result can be checked by numerical calculation. For simplicity, we consider the one spin-s 
system (n=l). A typical situation is shown in Fig.([T]). The agreement is spectacular. 

We can say a word on how the Bethe equations were solved. We first determine S(e) by solving the 
polynomial equation eq.© and then determine ip(z), eq. (fTT1) . by solving eq. (fT2l) . The Bethe roots 
are then obtained by solving the polynomial equation ip{z) = 0. 

The idea that the Bethe roots condense in the semi classical limit to form the cuts of the spectral 
curve goes back to [12]. It plays a very important role in the recent studies on the AdS/CFT 
correspondence in which it was greatly developed [THE]. Eq. (|26j) however seems to be new. 



7 Appendix: The XXX spin chain 



In the Gaudin model, the Bethe equations were shown to be equivalent to a Riccati equation eq. ([6]) . 
Moreover this equation itself determines the parameters S(ej), i.e. the eigenvalues of the commuting 
Hamiltonians. We show that this construction can be extended to the XXX spin chain. 



In the case of the XXX spin chain the Bethe equations take the form (see e.g. |15| ) 

and the corresponding generating function for the eigenvalues of the commuting Hamiltonians is 

m) = (* + ?)"n + (* - f ) " n (*) 

Let us introduce the polynomial 

M 

Q(\)= H(\-Vm) 

771=1 

Then the Bethe equations eq. (|27p can be rewritten as 
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Figure 1: Red dots are the Bethe roots ^ for the one spin system. Green dots are the branch 
points. The thin black curve is the solution of eq. (|26l) . (H = 1/30, s = 1/H, M = 4/H, highest 
energy state). ^3 



This means that the polynomial of degree N + M 



ih\ . , / . in 



A + yJ Q(X-iH) + \X- jj Q(X + ih) 

is divisible by Q(X). Hence there exists a polynomial t(X) of degree N such that 

■fj\N / ■f,\ N 

A + y) Q(A-a) + U-yJ Q(A + in)=t(A)Q(A) (29) 

This is Baxter's equation. The polynomial t(X) is the same as in eq. ()28p because that equation can 
be rewritten as 

hence the coefficients of this polynomial are the eigenvalues of the set of commuting Hamiltonians. 



Just as in the Gaudin model, it is interesting to introduce the Riccati version of eq. (|29|) . We set 

Q(A - ih) 



S(X) 

Then Baxter's equation becomes 



Q(X) 



(a + y ) N S(X) + (a - I ) N 5" 1 (A + ih) = t(X) 

This equation determines both S(X) and t(X). To find the equation for i(A), we expand around 
A = —ih/2 getting 

(e - ih) N S~ l (e + ih/2) = t(e - ih/2) - e N S(e - ih/2) 
Similarly, expanding around A = ih/2 we get 

(e + ih) N S(e + ih/2) = t(e + ih/2) - e N S' 1 (e + 3ih/2) 
Multiplying the two, we find 

ih\ , / ih\ 2 2 N 



t[e + y j t^e - yj = (^ + e'T + 0(e") (30) 

This is a system of iV equations for the TV + 1 coefficients of t(X) which determines it completely if 
we remember that t(X) = 2X N + 0(X N ~ 1 ). Eq. (|30p is the characteristic equation of the commuting 
Hamiltonians of the XXX spin chain. 

This construction can be generalized to the case of a spin-s chain. Baxter's equation reads 

(A + ihs) N Q(X - ih) + (A - ihs) N Q(X + ih) = t(X)Q(X) (31) 
and the Riccati equation becomes 

(A + ihs) N S(X) + (A - ihsf 5~ X (A + ih) = t(X) (32) 
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Taking s = 1 for instance, we expand around A = ih, A 



0, A 



ih getting 



(e + 2ih) N S(e + ih) 
(e + ih) N S{e) + (e - ih^S'^e + ih) 
(e-2^) JV 5- 1 (e) 



t( e + ih) + 0(e N ) 
t(c) 

t(e - ift) + 0(e*) 



from which we deduce 



i(e + ih)t(e)t(e - ih) = (e - i^T (e + 2ih) N t{e - ih) + (e + ih) N (e - 2ih) N t{e + ift) + 0(e JV ) 



Clearly for a spin-s, s > 0, the degree of the equation is 2s + 1. If however s < the equations 
generically do not lead to a finite degree equation as expected. 

In the semi classical limit H — ► 0, Hs —> s c i, eq. (132p tends to 



which is nothing but the spectral curve of the classical spin chain. 
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